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Abstract. While thermoelectric transport theory is well established and widely applied, it is not always 
clear in the literature whether the Seebeck coefficient, which is a measure of the strength of the mutual 
interaction between electric charge transport and heat transport, is to be related to the gradient of the 
system’s chemical potential or to the gradient of its electrochemical potential. The present article aims 
to clarify the thermodynamic definition of the thermoelectric coupling. First, we recall how the Seebeck 
coefficient is experimentally determined. We then turn to the analysis of the relationship between the 
thermoelectric power and the relevant potentials in the thermoelectric system: As the definitions of the 
chemical and electrochemical potentials are clarified, we show that, with a proper consideration of each 
potential, one may derive the Seebeck coefficient of a non-degenerate semiconductor without the need to 
introduce a contact potential as seen sometimes in the literature. Furthermore, we demonstrate that the 
phenomenological expression of the electrical current resulting from thermoelectric effects may be directly 
obtained from the drift-diffusion equation. 

PACS. 84.60.Rb Thermoelectric, electrogasdynamic and other direct energy conversion - 72.20.Pa Ther¬ 
moelectric and thermomagnetic effects 


1 Introduction 

Thermoelectricity is a mature yet still very active area of 
research covering various fields of physics, physical chem¬ 
istry, and engineering. The large interest in thermoelec¬ 
tric systems is mostly due to the promising applications 
in the field of electrical power production from waste heat 
as thermoelectric devices may be designed for specific pur¬ 
poses involving powers over a range spanning ten orders 
of magnitude: typically from microwatts to several kilo¬ 
watts. Further, thermoelectricity also provides model sys¬ 
tems that are extremely useful in the development of the¬ 
ories in irreversible thermodynamics nu. 

The discovery of the thermoelectric effect is usually 
attributed to Seebeck. In 1821, he published the results 
and analysis of his experiments aiming at establishing a 
magnetic polarization in a metallic circuit simply by per¬ 
turbing the thermal equilibrium across this latter [3]. More 
precisely, Seebeck described the appearance of a magnetic 
field within a closed electrical circuit made of two dis¬ 
tend offprint requests to: Y. Apertet 
a Email address: yann.apertet@gmail.com 


similar materials as the junctions between these materials 
were maintained at different temperatures. While Seebeck 
interpreted the observed phenomenon as a thermomag¬ 
netic effect, Oersted soon reexamined Seebeck’s work and 
showed that in this case the magnetic field was an indi¬ 
rect effect as it originated in the presence of an electro¬ 
motive force induced by the temperature difference [3], 
The proportionality coefficient between this electromotive 
force and the temperature difference across the system is 
the thermoelectric power, which has also been coined as 
“Seebeck coefficient”. 

The definition of the thermoelectric coupling has later 
been extended from that derived from the first experi¬ 
ments to both thermodynamic 0 and microscopic isimisi 

properties of materials. However, as of yet, there still 
is no clear consensus on its relationship with the vari¬ 
ous thermodynamic potentials and their variations (see, 
e.g., refs. amumnEniisi). Indeed as the terminol¬ 
ogy and conventions may vary from a discipline to an¬ 
other, say, e.g., solid-state physics and electrochemistry, it 
is not always straightforward to establish a clear distinc¬ 
tion or relevant associations between Fermi energy at zero 
or finite temperature, electrochemical potential, voltage, 
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Fig. 1. Determination of the Seebeck coefficient for a circuit 
composed of two dissimilar materials. 

Fermi level relative either to the conduction band mini¬ 
mum or to the vacuum, and chemical potential. 

In this article, we discuss the definition of the See¬ 
beck coefficient focusing particularly on the distinction 
between chemical and electrochemical potentials. First, in 
sect. [21 we address the experimental determination of the 
Seebeck coefficient in order to identify the quantities of in¬ 
terest. Next, the purpose of sect. [3]is to demonstrate that a 
clear physical picture of thermoelectric phenomena at the 
microscopic scale may be obtained on the condition that 
the potentials are carefully introduced. For this purpose, 
we review the standard definitions given in the literature 
to remove any confusion between the chemical and elec¬ 
trochemical potentials before we present and discuss our 
derivation of the Seebeck coefficient for a non-degenerate 
semiconductor. 

2 Experimental determination of the 
thermoelectric power 

The determination of the Seebeck coefficient traditionally 
involves components made of dissimilar materials, which 
we label A and B respectively. The two materials are com¬ 
bined to obtain two junctions as depicted in fig. [Q These 
junctions are then brought to different temperatures Tj 
and T 2 . An isothermal voltage measurement at a tem¬ 
perature T 3 , is performed between the free ends of the 
component B. The voltage thus measured is V 2 — V\ (this 
notation allows to clearly define a direction for the volt¬ 
age) and the Seebeck coefficient ciab associated with the 
global system, i.e. the couple AB, is defined as the pro¬ 
portionality coefficient between the resulting voltage and 
the applied temperature difference: 


The coefficient ctab, obtained for the whole circuit, is re¬ 
lated to the Seebeck coefficient of each material through 

m 

«AB = «B — CL A, (2) 

where «a and «b are the Seebeck coefficients of the ma¬ 
terials A and B respectively. 

From an experimental viewpoint, the presence of the 
material B A) is mandatory as it is associated with 
the probe’s wires (see, e.g., ref. [ill). However, if its See¬ 
beck coefficient ccb is sufficiently small to be neglected, 
the measurement may be used to determine directly the 
Seebeck coefficient of material A. In this case, one gets: 


Note the presence of a minus sign in the expression above: 
It is often overlooked in the literature but, fortunately, 
that omission is most of the time compensated by the ab¬ 
sence of a clear sign convention for the measured voltage. 

Let us now turn to the analysis of the measured quanti¬ 
ties. While the temperature is not subject to questioning, 
the voltage obtained from a voltmeter must be defined 
unambiguously. Indeed, it appears that its connection to 
the microscopic and thermodynamic properties of mate¬ 
rials has remained unclear for quite some time, leading 
Riess to publish in 1997, hence fairly recently, an article 
untitled “ What does a voltmeter measure ? ” m■ In that 
paper, Riess demonstrated that the voltage measured by 
a voltmeter between two points in a circuit is the differ¬ 
ence of electrochemical potentials jl at the two considered 
points divided by the elementary electric charge e, but not 
the difference between the electrostatic potentials p alone. 
The potential V might thus be defined as V = — jS/e. This 
result is recovered when one measures the voltage at the 
ends of a pn junction at equilibrium: While there is a built- 
in electric field associated with the depletion layer, the 
measured voltage remains zero. The Seebeck coefficient 
thus appears as a link between the applied temperature 
difference and the resulting difference of electrochemical 
potential between the two junctions. 

The simple technique presented here is not the only 
one used to determine the thermoelectric power of a given 
material. Indeed, since the measurement always involves 
a couple of materials, the absolute Seebeck coefficient of 
the second material has to be known accurately. To obtain 
this value, it is possible to use low temperature measure¬ 
ment to reach superconducting state where a = 0 and then 
derive higher temperatures values using the Thomson co¬ 
efficient that can be measured for a single material. For a 
detailed presentation of the Seebeck coefficient metrology, 
the reader may refer to the instructive review by Martin 
et al. [ 20 j . 

3 Relationship between the thermoelectric 
power and the electrochemical potential 


«AB 


V2-V1 

T2-Tff 


( 2 ) In order to better understand the influence of each po¬ 
tential, we identify the respective effects of temperature 
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bias, concentration difference, and electric charge, and we 
discuss the relationship between chemical potential, elec¬ 
trochemical potential and the band diagram of materials. 
We then derive the Seebeck coefficient in the simple case 
of a non-degenerate semiconductor to illustrate the con¬ 
tribution of each potential. 


3.1 Definition of the thermopower 

The Seebeck coefficient may be obtained from a micro¬ 
scopic analysis of the considered materials, with the local 
version of Eq. m, in open-circuit condition, i.e., with a 
vanishing electrical current: 


a = 


V/i 

eVT’ 


(4) 


where Jl and T are respectively the local electrochemical 
potential and temperature, defined at each point of the 
system. The notation V is associated with the gradient of 
each quantity. In the following, for the sake of simplicity, 
we consider a unidimensional system so that the spatial 
gradient reduces to its x-component: V, 7; . 


3.2 Distinction between the potentials 

Consider a semiconductor sample at thermal equilibrium 
and characterized by a spatially inhomogeneous doping. 
As the carrier concentration is nonuniform, a particle cur¬ 
rent takes place from the region of higher concentration to 
that of lower concentration: This is the diffusion process 
associated with the variation of the carriers’ chemical po¬ 
tential across the system. This type of electrical current 
is referred to as the diffusion current. The inhomogeneous 
electron population in the system thus generates an elec¬ 
tric potential difference and hence a built-in electric field 
which influences the electrons’ motion in such a fashion 
that it tends to curb the diffusion current. The electron 
motion driven by the built-in electric field is the drift cur¬ 
rent, which, at thermal equilibrium, exactly cancels the 
diffusion current, in accordance with the principle of Le 
Chatelier and Braun. In this case, the measured voltage 
across the system always remains zero and there is no 
net electrical current even if the system is short-circuited: 
The electric field associated with the electrical potential 
variation is obviously not an electromotive field. However, 
if the electrons are placed in a non-equilibrium situation 
caused by a thermal bias applied across the system, a non 
vanishing electric current may be obtained when the cir¬ 
cuit is closed. This current obviously stems from the un¬ 
compensated contributions of both the diffusion and drift 
of charge carriers, and it is traditionally related to the 
gradient of the temperature and to the gradient of the 
electrochemical potential. 

The electrochemical potential Jl of a population of elec¬ 
trically charged particles is the sum of a chemical con¬ 
tribution /z, the chemical potential, and of an electrical 
contribution H3: 



Fig. 2. Energy levels in an n-type semiconductor highlighting 
the notations used in this article (adapted from ref. [8]). The 
energy Eg refers to the bandgap energy. 


M = At + Me- (5) 

Note that the quantities we just referred to as potentials 
are actually energies. The electrical contribution /i e may 
be expressed as a function of the electrostatic potential 
p (a genuine potential contrary to Jl and p) so that the 
electrochemical potential reads: 

Jl = p + qp, (6) 

where q is the electrical charge of the considered particle. 
When used in solid state physics, these quantities have 
to be related to an energy band diagram. This correspon¬ 
dence may be found for example in the book of Heikes and 
Ure [5]: Considering the example of an n-doped semicon¬ 
ductor, the electrochemical potential Jl corresponds to the 
Fermi level, the electrostatic energy —ep corresponds to 
the energy level of the bottom of the conduction band 
while the chemical potential p corresponds to the dif¬ 
ference between these two quantities and is often called 
Fermi energy These notations are summarized on fig. ©. 
The difference between Fermi level (p) and Fermi energy 
(p) was already highlighted by Wood [3]: 11 The differ¬ 
ence between the Fermi energy and the Fermi level should 
be noted. The Fermi energy is generally measured from 
the adjacent conducting band edge (valence or conduction 
band for holes or electrons, respectively), i.e. a reference 
level which may vary in energy, whereas the Fermi level is 
measured from some arbitrary fixed energy level ”. This last 
remark stresses the importance of the choice of an energy 
reference, which is a key parameter: To express energies 
in a semiconductor, the bottom of the conduction band 
is often used as the reference uni; however, for studies 
of non-equilibrium phenomena such as thermoelectricity, 
it is mandatory to define an arbitrary fixed energy refer¬ 
ence independent of the position within the material since 
both Jl and p may vary along the system. It seems the 
only way to correctly describe the relative displacement 
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T 3 Ti T 2 T 3 

Fig. 3. Schematic illustration (adapted from Ref. [21]) of the 
variations of the bottom of the conduction band, Ec, the top 
of the valence band, Ey, the Fermi level, Ef, and the vac¬ 
uum level just outside the material, es, all along the circuit 
depicted in Fig. [l] The slopes of the lines have been greatly 
exaggerated for clarity, and band bending at the interfaces has 
been neglected. 


of these energies. Note that the vacuum level infinitely far 
from the system, Eoo, might be a good and meaningful 
energy reference. 

Figure [3] illustrates the variations of the different en¬ 
ergies around the circuit depicted in Fig. [T] in the case of 
semiconductor materials. It highlights the difference be¬ 
tween the slope of the bottom of the conduction band and 
the slope of the Fermi level: The variation of the chemical 
potential thus differs from the variation of the electro¬ 
chemical potential. Distinguishing these two energies is, 
therefore, crucial to properly evaluate the Seebeck coeffi¬ 
cient. Figure [3] also displays the vacuum level es just out¬ 
side the material (different from .Eoo)- This vacuum level 
is related to the bottom of the conduction band through 
the affinity \ of the material. The discontinuities in eg 
at the interfaces might be seen as contact potentials. On 
the contrary, the Fermi level Ep is continuous along the 
system, even at the interfaces. Its variation however un¬ 
dergoes a sudden change at the interface, reflecting both 
changes in temperature gradient (assumed constant in a 
given material) and in Seebeck coefficient from a material 
to an other. The thermopower is indeed associated with 
bulk material but not to interfaces. A similar figure for a 
system made of metals can be found in Ref. [21 1 . 


3.3 From potentials to thermoelectric power: the 
illustrative case of a non-degenerate semiconductor 


„(T) = JVexp = "exp (AY) , (7) 


with 


^ 27rTO e fffc B r ^ 3/2 


( 8 ) 


and where Ec is the energy level of the bottom of the 
conduction band, m e g is the electron effective mass, fc B 
is the Boltzmann constant and h is the Planck constant. 
The Seebeck coefficient is associated with non-equilibrium 
phenomena, and, as such, it is tightly linked to transport 
properties of electrons inside the material. To take account 
of these properties, we build on the drift-diffusion equation 
used to obtain the net electrical current density J x : 


J x — cnM n £ x T cD n ^ x n, (9) 

where M n and D n are the electron mobility and diffusivity, 
and where the electric field £ x is related to the energy level 
Ec through: 


e x = = YjeE .£. ( 10 ) 

q e 

At first, we assume a situation where the electron diffusiv¬ 
ity D n does not depend on the other parameters, including 
the position. The variation of D n will be discussed further 
below. 

The Seebeck coefficient is obtained setting J x = 0. 
However this current density should be related first to 
S7 X T and V X J1 rather than to £ x and V x n. To do so, we 
evaluate the gradient of the electron density given by eq. 
0 considering that Eq , ^ and T may vary along the mate¬ 
rial. This approach is seldom found in the literature as one 
often sets Eq = 0 , thus considering the bottom of the con¬ 
duction as the reference everywhere in the nonequilibrium 
system. As already stressed, this viewpoint is misleading 
for thermoelectric phenomena. From eq. 0 , the gradient 
of electron density reads: 

Q 77 \7 r f 1 77 

S7 x n=-^^ + — [T{S7 x Ji- e £ x )-^ x T]. ( 11 ) 

We then use this equality along with Einstein’s relation 
between the electron mobility M n to the electron diffusiv¬ 
ity D n : 


M n __ e 
D n /c b T 

to modify eq. m as follows: 


( 12 ) 


We emphasise the importance of the distinction between 
ji and /i on the derivation of the thermoelectric power us¬ 
ing the example of a non-degenerate semiconductor doped 
with electrons. In this case, the expression of the carrier 
concentration n is rather simple: 


Jx — CTlAJ n 

Now, setting J x 
eq. 0, we find: 




3 

2 ~~ ksT 


V X T 


(13) 


ising the definition given 


m 
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Scattering mechanism 

Exponent s 

Acoustic phonon 

-1/2 

Ionized impurity (strongly screened) 

-1/2 

Neutral impurity 

0 

Piezoelectric 

+ 1/2 

Ionized impurity (weakly screened) 

+3/2 


Table 1. Values of the exponent s for different scattering 
mechanisms (adapted from Ref. 1221 1. 


a = — 


k B 


e 


3 

2 


k B T 


(14) 


with a constant electron diffusivity, which is the expected 
expression for a non-degenerate semiconductor. Further, 
this result may also be interpreted by looking at the net 
thermal energy transported by each carrier transported 
inside the material, i.e., qll , where 77 is the Peltier coef¬ 
ficient m ■ For electrons, this energy is 


- e77 = -jk B T - (J,, 


(15) 


since it corresponds to the energy above the Fermi level 
J1 and hence to the sum of the average thermal energy for 
free electrons and of the energy between the Fermi level 
and the bottom of the conduction band, i.e., — \i. We thus 
recover the second Kelvin relation relating the Seebeck 
and Peltier coefficients: 77 = aT. 


in the conduction band and may thus be approximated by 
its average value, i.e., 3/2k B T. Replacing r in eq. (fl7)l . we 
obtain: 


V x (£>„) = ^(l + s)V x (T). (18) 

Finally, inserting eq. (TTH1) and eq. m in eq. © yields: 


J x = enM n 


k B 

e 


5 n 

2 +S ~k^T 



and consequently: 


cr = — 


k B 


e 


5 fJL 

2 +S ~k^T 


(19) 


( 20 ) 


The contribution of the diffusivity gradient to the thermo¬ 
electric power is — (1 + s)k B /e and hence depends only on 
the scattering parameter s. This term has also been recov¬ 
ered by Cai and Mahan m using a Boltzmann equation. 
Note that this term was also introduced by Ioffe [23] with 
the notation old- However, Ioffe used a different power 
law: He assumed that the carrier’s mean free path l is 
proportional to (E — EqY■ He consequently found that 
old = —(1/2 + r)k B /e. This discrepancy is quite easy to 
understand since r is proportional to l/y/E — Eq. 


4 Discussion 


3.4 Taking into account diffusivity variation 

If we relax the assumption of constant diffusivity D n , this 
latter becomes a function of the spatial coordinate x and 
we end up with the so-called Stratton equation [22] : 

Jx = cnAJ n S x + x (D n n) 

= enM n £ x + eD n X7 x ( n ) + ensJ x (D n ) , (16) 

It corresponds to a more general form of the drift-diffusion 
equation, which contains a third contribution to the car¬ 
rier motion, directly linked to the gradient of diffusivity 
along the system. To evaluate its effect on the thermo¬ 
electric power, we may reexpress it as a function of the 
temperature gradient using the relation between the dif¬ 
fusivity D n and the relaxation time of the carriers r. Since 
M n = eT/m e g , the Einstein relation reads: 

D n = k B T -. (17) 

mefi 

To keep the calculations on an analytical level, we assume 
that we deal with low-energy conduction electrons, and 
we express the relaxation time using a power law of the 
form: r oc (E — Ec) s , where E is the total energy of the 
carrier and s is a characteristic exponent depending on 
the scattering mechanisms [22]. Some typical values for 
this exponent are given in Table [Q Note that the energy 
E — Ec corresponds to the thermal energy of the carriers 


4.1 An unusual derivation 

While the result given in Eq. m is well-known, its deriva¬ 
tion presented here is quite original. Indeed, it was di¬ 
rectly obtained from the drift-diffusion equation. Thus, 
the phenomenological equation associated with thermo¬ 
electric transport: 

J x = a — - raVJ, (21) 

e 

where cr = enM n is the electrical conductivity, is identi¬ 
cal to eq. m (or to eq. (fl3l) depending on the hypothesis 
made). This latter appears as a modified form of the drift- 
diffusion equation, which accounts for the couple of vari¬ 
ables [ju, T], or more precisely their gradients, rather than 
the traditional couple [n, if]. This modification puts forth 
the fact that the first term of the right hand side of eq. 
m does not correspond any longer to the genuine local 
form of Ohm’s law since it does not involve the electrical 
field £ x . In this case, the true electromotive force is given 
by the gradient of the electrochemical potential as carriers 
experience both diffusion and effects of the electric field. 

The simple derivation of Eq. m has been allowed by 
the use as a reference of a fixed energy level, arbitrarily 
chosen but independent of the position along the material, 
rather than the bottom of the conduction band. This ap¬ 
proach demonstrates that a particular knowledge of both 
V x n and £ x is not mandatory to determine the Seebeck 
coefficient of a non-degenerate semiconductor. Indeed, in 
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this case, these two contributions to the electrochemical 
gradient seem to always compensate each other in such 
a way that the resulting electromotive power is indepen¬ 
dent of specific assumptions, for example a constraint on 
the carrier concentration. Equation (1191) is thus valid for a 
wide range of temperatures: It is correct for the extrinsic 
regime, i.e., when the carrier concentration is fixed by the 
concentration of impurities, but it also remains valid in the 
freeze-out regime and in the intrinsic regime where addi¬ 
tional carriers are thermally generated. However, in this 
latter regime, the electron hole contribution to thermo¬ 
electric power should also be considered as these minority 
carriers may no longer be negligible. One may also refer 
to Ref. HU in which the authors derive the Seebeck co¬ 
efficient focusing only on potentials and the electric field 
rather than using a statistical approach. 

4.2 Link with non-equilibrium thermodynamics 

While Eq. ETl) is widely used in solid-state physics, its 
formulation is slightly different in non-equilibrium ther¬ 
modynamics as general forces are traditionally computed 
from the gradients ^\7 x jl and V x (1 /T), instead of V x jl 
and V X T HZH23- So, one should then rewrite Eq. m to 
get: 

J x = ^ V x [i + craT 2 ■ V x ^ . ( 22 ) 

With this form, it is possible to identify each term with 
the canonical expression m , 

— Jn = —^ = Ln ■ —V x n + L 12 ■ V x > (23) 

to recover the expressions of the kinetic coefficients in the 
thermoelectric case, i.e., in = oT/e 2 and L\ 2 = aaT 2 /e. 


i.e., to eq. ED- This latter term is identical to the so-called 
Kelvin formula of the thermopower [13]. As discussed by 
Shastry HU, this contribution “captures the many body 
density of states enhancements, while missing velocity and 
relaxation contributions ' 1 '’. It thus justifies the introduction 
of the coefficient S to take into account dynamical effects. 
We believe however that this coefficient should not be pre¬ 
sented as effective since it does not reflect the appearance 
of the electromotive force due to the temperature gradi¬ 
ent. It represents only one of the possible contributions 
to this electromotive force. From a practical viewpoint, it 
has recently been demonstrated that the contribution to 
the thermopower from the diffusivity gradient might be 
significant m- 

4.4 On the contact potentials 

Finally, we want to point out the inappropriate use of the 
contact potentials in the derivation of the thermoelectric 
power sometimes found in the literature. For example, in 
Ref. HSi Ioffe obtains Eq. (flTJl) splitting the Seebeck co¬ 
efficient into three separate terms, one being ap> while 
the two others, a n and a are associated respectively 
to concentration gradient and to the “temperature depen¬ 
dence of the contact potential”. However, as demonstrated 
later by Chambers HU: contact potentials are irrelevant 
to thermoelectric effects. This latter term is indeed intro¬ 
duced only to compensate the erroneous expression of a n 
stemming from the confusion between ip and p/e. From an 
experimental viewpoint, contact potentials are irrelevant 
since they cannot be probed by a voltmeter: As depicted in 
Fig. El these energy discontinuities concern only the bot¬ 
tom of the conduction band Eq (or identically the vacuum 
level just outside the material eg) but not the Fermi level 
Ep. 


5 Conclusion 


4.3 On the so-called effective Seebeck coefficient 


Let us now turn to the previous analysis of the thermo¬ 
electric power in non-degenerate semiconductor. In ref. 
m, Mahan introduces an effective Seebeck coefficient S, 
distinct from the genuine thermoelectric power eq. (El) ob¬ 
tained from measurements. In a subsequent article with 
Cai [12] . this effective coefficient is presented as the ra¬ 
tio between the electric field and the temperature gradi¬ 
ent. These two different Seebeck coefficients are related 
through the following relation m : 


= S 


^ 1 


e 



(24) 


A comparison of eq. (1161) with the equation (20) of ref. 
m leads to identify the effective Seebeck coefficient to 
the contribution of the diffusivity gradient, i.e., ap>. The 
second term of the right hand side of eq. (12H) should then 
be associated with the assumption of constant diffusivity, 


In this article, we have discussed the definition of the 
thermoelectric power with a special emphasis on its re¬ 
lationship to the electrochemical potential. A proper con¬ 
sideration of all potentials inside the material has led to 
demonstrate that the phenomenological equation for the 
electrical current involving thermoelectric coefficients may 
be derived directly from the drift-diffusion equation. We 
also shed light on the physical interpretation of the effec¬ 
tive Seebeck coefficient defined by Mahan, showing that it 
is actually related to the gradient of diffusivity along the 
system. 
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